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, We demonstrate that in the strong coupling limit (the superconducting gap 

A is as large as the chemical potential /i), which is relevant to the high-T c 
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I/-) ■ superconductivity, the correlation corrections to the gap and critical temperature 

l> . 

are about 10% of the corresponding mean field approximation values. For the 

m 

weak coupling (A -C /i) the correlation corrections are very large: of the order 
of 100% of the corresponding mean field values. 
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I. INTRODUCTION 

In recent papers we demonstrated the very strong d-wave pairing between dressed 
quasiholes in the t — J model induced by spin-wave exchange. The pairing gives the critical 
temperature in a reasonable agreement with experimental data for Cooper Oxide Supercon- 
ductors. For calculations we used BCS-like mean field approximation for dressed quasiholes. 
A similar approach to the pairing of dressed quasiholes has been used by Dagotto, Nazarenko, 
and Moreo 0. The important difference is that in the hole- hole interaction was derived 
from the parameters of the t-J model, while in || it was introduced ad hoc with magnitude 
adjusted to fit experimental data. The typical value of the gap A obtained in the papers Jl] [| 
is of the order of the chemical potential fi. This strong pairing inspires a natural question : 
how strongly the hole-hole correlations influence upon the mean field result? The purpose of 
the present work is to investigate this problem. We consider the simplified model instead of the 
sophisticated t-J one. The model under consideration describes the two-dimensional fermions 



with quadratic dispersion and contact attraction. This model inherits the main, relevant to the 
problem, property of the t-J model: it permits to consider strong pairing. Therefore one can 
suppose the obtained result to be rather general. We study the dependence of correlations on 
the intensity of attractive interaction. Our conclusion is that for strong coupling A > /x the 
correlation corrections are not large: about 10% of the corresponding mean field value. This 
is a very surprising result since one could expect that strong interaction causing the pairing 
might make the correlation correction to be very large. Our conclusion is in agreement with 
the one which was recently obtained in finite-cluster numerical study of t — J model by Ohta, 
Shimozato, Eder, and Maekawa ||. 

The other surprising result of the present work is that in the weak coupling limit A <C /j 
the correlation corrections are very large: the renormalizations of mean field values are about 
100%. This is due to the specific ultraviolet behaviour in the two dimensional theory with 
attraction. 

We are mainly interested in the regime A > /i when the system is close to a smooth crossover 
from a state with large, overlapping Cooper pairs to a Bose condensate of composite bosons. In 
the mean field approximation this crossover has been studied in three dimensions by Legget || 
and by Nozieres and Schmitt-Rink ||. A similar problem for two dimensions reveals interesting 
features considered by Randeria, Duan, and Shieh J/J. In the present work we investigate the 
correlation corrections to the results presented in Ref. 0] 



II. MEAN FIELD APPROXIMATION 

Consider the Hamiltonian of fermions with spin 1/2 and contact attractive spin independent 
interaction 

H = J2 2m a ^ ak ' a ~v ^2 4 1+ k 2 -k3,T a k 3 a ak 2a a ki,T> (1) 

k,cr ki,k 2 ,k 3 

where summation over k is restricted in two-dimensional plane, V is the area of the plane, a = 
±1/2 =|, | is a projection of the usual 3-dimensional spin 1/2. Consider first the Schroedinger 
equation for a two-particle bound state with zero total momentum 

Xk = TT ~^ K~/ • (2) 

V ~ E a - p 2 /m 

The sum here is logarithmically ultraviolet divergent, and therefore one has to introduce the 
ultraviolet cutoff E\ = A 2 /2m. For the t — J model the parameter A is of the order of inverse 
lattice spacing. Solution of Eq.(Q) is straightforward 



Xp = const, E a = -2E A e~^ m9 . (3) 

Now let us consider a many body problem with fixed number density of particles 5. The Fermi 
energy is: Ep = irS/m. The BCS equation for the pairing at fixed chemical potential fi is of 
the form 

A k = ^-y^Ptanh^, (4) 



where e p = ^Ji] p + A p , and r/ p = p 2 /2m — /x. Similarly to the Schroedinger equation we have 
to introduce the ultraviolet cutoff A. Solution of the BCS equation for zero temperature is 
straightforward. Assuming that E F <C E^ one gets 



A = pE F \E a \ = 2^E F E A e~ 2 ^ m9 , (5) 

1 A 2 ' 



/ 1 A z \ 

fi = E F -\E a \/2 = E F (l--g r ) 



F ' 

where E a is the binding energy (|3|) of the two particle bound state. This solution obtained by 
Randeria, Duan, and Shieh H] gives a smooth crossover from the BCS limit (// ~ Ep) to the 
Bose condensate of composite bosons (// < 0). In the present work we concentrate on the case 
of positive \i because, in our opinion, it is relevant to the realistic high-T c superconductors. 

The critical temperature T c can be easily found from Eq.([|) if we define it as a point where 
the gap vanishes. In the weak coupling limit (T c Ep) we have the usual BCS relation: 
A(T = 0)/T c = 1.76. Numerical solution of Eq.(^j) shows that even in the strong coupling limit 
this ratio remains very close to the BCS value. For example at A(T = 0)//i(T = 0) = 3 which is 
equivalent to A(T = 0)/E F m 1.5 one finds A(T = 0)/T c = 1.68. We assume that the density 
of particles is fixed. One has to remember that under this condition the chemical potential 
is a function of temperature: fi(T = 0) 7^ /i(T = T c ). For strong coupling the temperature 
dependence of the chemical potential is not negligible. 

Equation (|J) as well as solution (|5]) describe the mean field approximation. Now let us 
consider correlation corrections. We use the conventional Gorkov-Nambu technique, see e.g. 
book M. Consider first the case of zero temperature. 

III. CORRELATION CORRECTION TO THE GAP AT ZERO TEMPERATURE 



The normal G(l,2) = -i(T[^(l)^t(2)]) and the anomalous 2^(1,2) = -i(T[^ t (l)^ t (2)]) 
Green functions obey the usual Dyson equations || 



G(p) = [1 + G(jp)± xl (jp) + F^p)t 20 (p)}G (p) (6) 
F\p) = [F^p)t u (-p) + &p)t 02 (p)]G (-p), 

where p = (e, p), Gq(p) = [e — r] p + iO ■ sign^-p)]^ 1 . In the first order of perturbation theory 
the normal self-energy operator E^ is given by the diagrams presented in Fig.l. With the 
interaction ([!]) the self-energy is momentum independent, and therefore it gives only a 
correction to the chemical potential /x. The first order anomalous self-energy operator E^ is 
given by the diagram Fig. 2. It is equivalent to the BCS mean field approximation. Solution 
of Eq.(||) with E 2 o = E20 is of tne form (see [§): G a/3 (p) = 5 a pG{p), F^p) = g a pF\p), 
E(20)a/3(j>) = ga/3^20 {p) , where 6 a /3 and g a p are standard symmetric and antisymmetric spin 
matrices, and 



e p + iO e + e p - iO 



Ft(e, p) = -u p v p ( 1 —— - 1 ) (7) 

E (1) - A 
^20 _ ^ 



with Up, Dp = ± r/ p /e p ). The gap A is given by (|). 

(2) 

In the next order of perturbation theory the normal self-energy operator E^ (p) is repre- 
sented by the diagrams in Fig. 3. Note that that inside E^ 2 ) we use not the bare Green functions, 
but the "dressed" ones. They take into account the self-energy corrections in accordance with 
Eq.(|5]). This approach is widely used in many-body problems when a correlation correction can 
be significant, see for example Ref. ||. The self-energy corresponding to diagrams in Fig. 3 can 
be easily evaluated. 

4? (e, P) = (jX E [vlul - (um) (u 2 v 2 )} f Ul 



y) k f^ 2 Ve-ei-e 2 -e3 -e-ei-e 2 -e 3 . 

The summation here is carried out over and k 2 . The momentum k 3 is defined as 
k 3 = p + k x + k 2 . The second order anomalous self-energy operator is given by the diagrams 
presented in Fig. 4. The calculation gives 

E$(e,p) = (f) 2 E Usvsivful-Miu^f + ) (9) 

VV ' ktt 2 \e-e 1 -e 2 -e 3 -e- e l - e 2 - e z J 

Similar to (|8|) the momentum k 3 in this sum is equal to k 3 = p + kx + k 2 . 

Note that calculating the first order anomalous self-energy operator Fig. 2 we also have to 
use the exact Green function. Hence this operator is proportional to 



where is the "dressed" anomalous Green function, not just a first order one. 
Equations @ can be rewritten as 

[G \p) - X u (p))G(p) + (gE + Eg)Ft(p) = 1 (11) 
[G \-p) - E u (-p)]Ft(p) - (</H + Eg)G(p) = 0. 

We recall that is momentum independent and therefore is completely absorbed into the 
chemical potential. The second order self-energy (|8|) is logarithmically divergent at large k 2 . 
However, the diverging part is independent of e and p, and therefore can also be absorbed into 
the chemical potential. This freedom permits us to renormalize the self-energy imposing the 



condition £n(e = 0, |p| = p M ) = at p^ = v2m/I. We can use also the usual linear expansion 
near the point e = 0, |p| = p^: £n(e, p) = ^f^e + ^f^Vpi assuming that e ~ i] p ~ A. This 
expansion is certainly valid for the weak coupling A -C Ep. We have verified numerically using 
(|8|) that the expansion remains valid with a reasonable accuracy for the strong coupling A ~ Ep 
as well. Now we can easily find the solution of the Eq . (|Il"l) 

G{e,j>)nz( 5 + " l 



e — e p + iO e + e p — iO 

Ft(e, P ) » -Zu p vJ 1— - } ) (12) 

Ve — e p + 2(J e + e p — «0/ 



where Z = [1 — ^f 11 ] -1 is the quasiparticle residue, and the renormalized dispersion is e p 



?7p + A^ with 

A p = z(gE + X®(e,p)) (13) 



Numerical computations show that the dependence of A p on e and p at e ~ Ep and p ~ p^ is 
actually rather weak. But anyway, the gap depends on energy and momentum and therefore 
it is convenient to introduce A = A e= o iP=p , and to represent the second order anomalous 
self-energy as S^q (e, p) = —A ■ <r(e, p). From Eqs.@ and @ one finds 

f— ) =-(^) 2 R 1 (A/ l i) 
\ de J e=o,p=p„ V 27r / 

^ = r^ViiCA/rt, (M) 

2 



^ P / e =o, P =p M V 2tt / 



a(e = 0,p = p,) = (^-) Rs(A/f,). 



The functions Ri, i = 1,2,3 depend only on the ratio A//i. One can easily find that in the 
weak coupling limit (A <C ji): i?i(0) = const, #2(0) = const, i? 3 (A//i) In (///A). Results of 
numerical computations of at arbitrary A//x are presented in Table I. They show that the 
corrections due to En are negligible not only in the weak coupling limit when gm/2ir <C 1, 
A//x <C 1, but remain small for the strong coupling, gm/2^ ~ 1, A//i>l, as well. Therefore 
we can neglect En and consider only the anomalous self-energy operator a. The smallness of 
En corrections results in the fact that Z « 1. Therefore the exact Green functions (0) have 
the form similar to the Green functions (0) in the mean-field approximation. We used this fact 
when evaluated Eqs.(^),(^). 

In order to find a relation between the gap A and the coupling constant g we have to 
substitute the solution (p~2|) , (p~3|) into the self-consistency condition fllCf). Then we get 



V^2e p V^J 27rz(e-e p +z0)(e + e p -z0)' 1 ; 

The last term here gives the correction to the BCS mean field equation (f|). Assuming that the 
correction to the mean field value of the gap A m f is small (5 A = A — A m / <C A) we find from 



(ED 



6A 



_fgm\ 
A \2n) 

where L depends only on the ratio A/u 



2 



L(A//i), (16) 



L(A/a) = - , 2 * V / — ^ -. (17) 

1 ' ' \gm) i + Ji + A*/ia*V^J 27rz(e-e p + i0)(e + e p -i0) 1 ; 



Numerical computation of L is straightforward. It is convenient to integrate in (|17|) along the 
imaginary e axis. Results are presented in Table I. In the weak coupling limit (A <C fi) L can 
be easily calculated analytically with logarithmic accuracy: L « ln 2 (/i/A). Using flT|) and (|) 
we find the gap in this limit 



A « A 



m/ 



^In^'l-A fi f ln(^/A) ~ ! 



V2tt A 



A 



mf 



1 



Vln(2 v / ^El/A) 



;is) 



Thus the correction to the mean field value is very large. The reason for the large correction 
is simple. The gap is proportional to A oc e~ 2n / m9 . Practically we have calculated the renor- 
malization of the coupling constant g. A small correction to g gives a large correction to A 
when exponent 2ir/mg is large. Certainly in this situation the third order self-energy can give 
substantial contribution as well. 



For application to high-T c superconductivity we are more interested in the strong coupling 
limit A > il. Surprisingly in this case the correction to the mean field approximation is small. 
For illustration let us set E^/Ep = 20. With this ratio fixed one can easily find from (H) the 
value of gm/2iT as a function of A//i. After substitution of this value into with L from 
Table I we find <5A/A. Results of these calculations are presented in Table II. In our opinion 
the values of parameters A//i ~l-3 and E^/Ep ~ 20 correspond qualitatively to the t — J 
model describing high-T c superconductors We see from Table II that in this region the 

correlation correction SA/A is about -10%. 

We conclude that the BCS-like mean field approximation for the pairing of dressed holes in 
the t-J model is justified with the accuracy |<5A|/A ~ 0.1. It is worth to note that the situation 
when correlation corrections to the Hartree-Fock approximation are more important for weak 
coupling than for strong coupling is well known for a number of many-body problems in nuclear 
and atomic physics. 



IV. THE CRITICAL TEMPERATURE 

The correction to the critical temperature may be found in a way similar to the above 
developed approach for the correction to the gap at zero temperature. We have to solve Eqs. (§) 
and find a point where the gap vanishes. For finite temperature the energy in these equations 
is equal to e = i(, s , where £ s = 7rT(2s + 1), s = 0, ±1, ±2, ... is the Matsubara frequency. 
Integration over energy inside any loop should be replaced by summation over Matsubara 
frequencies. The mean field approximation is equivalent to the account of the diagram Fig. 2 for 
the self-energy. The solution of Eqs.(|5]) is of the form: G a /3(p) = 8a/3G(p), Fapip) = —Qa^ip), 
S(2 )a/3(p) = g a ^2oip), where 

G(i£s,p) = t— , 

*?* - V P 

Ss ~r '/p 
^20 — ZA - 

We neglect here all powers of the gap higher than one. Self-consistency condition for the 
diagram Fig. 2 gives the critical temperature in the mean field approximation. 

The second order normal self-energy operator is given by the diagrams presented in Fig.3a,c. 
The contributions of the diagrams Fig.3b,d,e,f are proportional to A 2 and therefore can be 
neglected. The calculation gives 



sS?(^,p) = E Hvi) ~ n{r]2)Mvi - m) - n(m] . (20) 

Here 71(77) = [1 + exp^/T)] -1 is a Fermi-Dirac function. The summation is carried out over 
k x and k 2 . The momentum k 3 is equal to k 3 = p + k x + k 2 . Similar to the case of zero 
temperature the detailed analysis demonstrates that the normal self-energy operator is small 
in the strong coupling limit as well as in the weak coupling one. Therefore below we neglect 
£11. 

The second order anomalous self-energy operator is given by the diagrams Fig.4a,c,d,e. The 
contribution Fig.4b is proportional to A 3 and therefore is neglected. The calculation gives 

V (2) 
^20 



l2) (7e,p) = -Aa T (^,p) 



a at ry) = (£_Y v K^) - n ^)\ ( fcM - n fa - V2)] _ [n(-m) - n(m - (21) 

TKSsiP) \v) £t 2r/ 3 I ^ s + 771-772-773 76 + ^71-772 + 773 P 1 ' 



where, as above, k 3 = p + ki + k 2 . After substitution of £ 20 into @ we find 

c + % 

where E* = T/FEsEp^siP) is given by the diagram Fig.2. Similar to (|T5|) this gives the 
equation for the critical temperature 

1 = L V J_ tanh l-lr c Vy ^ f • (23) 

The last term here is the correction to the mean field equation (f|). Assuming that the correction 
to the critical temperature is small we find from (E3|) 

f = (24) 

The function Lt depends only on the ratio T c //i c , where \i c is the chemical potential at the 
critical point, and 

mtm = (ii)' ' , 9T S T ^ a -wr4- < 25 > 

\gmj l + tanh(/i c /2T c ) V ^7 H + % 
Numerical computation of Lt is straightforward. Results are presented in the last column of 
Table I. We present Lt as a function of A(0)//x(0) (the gap at zero temperature over the chem- 
ical potential at zero temperature). It is possible to do so because T c j fi c itself is a function of 
A(0)///(0). In the weak coupling limit Lt can be easily calculated analytically with logarithmic 
accuracy: L T w ln 2 (/z c /T c ) « ln 2 (i?^/A(0)). From (|24]) and @ we find the critical temperature 
in this limit 



- c(m/) 



i ( \n(Ep/A) 



(26) 



Vln(2/E^EI/A), 

where T c ^ m ^ is the critical temperature in mean field approximation, and A = A(0) is the 
gap at zero temperature. Comparing fl2"E| ) with (|T^) we see that A and T c have the same 
renormalization factor. Therefore the BCS relation A/T c « 1.76 is preserved despite of the fact 
that the renormalizations of the mean field values are about 100%. 

Comparing L and Lt from Table I we see that in the strong coupling limit (A > /i) 
the correction to the critical temperature is larger than the correction to the gap at zero 
temperature. Nevertheless the correction remains small. Consider the same example as for 
zero temperature: E^/E F = 20. The Table II gives 5T C /T C as a function of A(0)//z(0). We see 
that at A(0)//i(0) ~l-3 the correlation correction ST C /T C is about -15%. 



V. CONCLUSION 



We consider the correlation corrections to the BCS mean field pairing in the two dimensional 
case. It is found that for the strong pairing (A > fi) the correlation corrections are not large: 
about 10% of the corresponding mean field value for the set of parameters relevant to t — J 
model describing high-T c superconductors. The small values of the corrections is explained 
qualitatively by the fact that the energy of virtual excitations becomes higher with increase of 
the pairing. We conclude that the BCS mean field approximation is reasonably justified for 
description of the dressed quasiholes pairing in the t — J model. 

Surprisingly for the weak coupling limit (A <C /x) the correlation corrections are very large: 
the renormalization of the mean field values is about 100%. The large correction results from 
the exponential dependence of the superconducting gap on the coupling constant which makes 
a small correction to g to give a significant contribution for A. 
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TABLES 



TABLE I. Numerical values of dimensionless functions Ri, L, Lt for different values of 
A(T = 0)//i(T = 0). We present also the corresponding values of A(T = 0)/Ep. The functions 
Ri,B,2 give normal self-energy operator En, see Eq.(14). Their small values permit one to neglect 
corrections caused by En and consider only corrections caused by £20- The functions L and Lt 
describe the correlation corrections to the gap and to the critical temperature, see Eqs.(16),(17) and 
(24),(25). 





A/Ep 


Ri 


R2 


R3 


L 


Lt 


3 


1.44 


3.4- 10~ 2 


1.9 • 10~ 2 


0.11 


0.27 


0.46 


2 


1.23 


5.0 • 10~ 2 


2.9 • 10~ 2 


0.15 


0.35 


0.59 


1 


0.83 


9.8 • 10~ 2 


5.8- 10- 2 


0.31 


0.65 


0.98 


0.5 


0.47 


0.19 


0.10 


0.66 


1.4 


1.9 


0.1 


0.1 


0.46 


0.17 


2.1 


6.6 


7.2 


0.01 


0. 01 


0.66 


0.24 


4.5 


23 


23 



TABLE II. The correlation corrections to the gap at zero temperature and to the critical temper- 
ature for different values of A(T = 0)//i(T = 0). The ultraviolet cutoff is fixed: E A /E F = 20. 



A//x 


gm/2-K 


-5A/A 


ST C /T C 


3 


0.55 


0.08 


0.14 


2 


0.51 


0.09 


0.15 


1 


0.42 


0.12 


0.18 


0.5 


0.34 


0.16 


0.22 


0.1 


0.22 


0.33 


0.36 


0.01 


0.15 


0.50 


0.50 



FIGURE CAPTIONS 

FIG. 1. Normal self energy operator in the first order of perturbation theory E n \ 

FIG. 2. Anomalous self energy operator in the first order of perturbation theory £20 • 

(2) 

FIG. 3. Normal self energy operator in the second order of perturbation theory E n . 
FIG. 4. Anomalous self energy operator in the second order of perturbation theory £ 
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